home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / spotf2.z / spotf2
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. SSSSPPPPOOOOTTTTFFFF2222((((3333FFFF))))                                                          SSSSPPPPOOOOTTTTFFFF2222((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SPOTF2 - compute the Cholesky factorization of a real symmetric positive
  10.      definite matrix A
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE SPOTF2( UPLO, N, A, LDA, INFO )
  14.  
  15.          CHARACTER      UPLO
  16.  
  17.          INTEGER        INFO, LDA, N
  18.  
  19.          REAL           A( LDA, * )
  20.  
  21. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  22.      SPOTF2 computes the Cholesky factorization of a real symmetric positive
  23.      definite matrix A.
  24.  
  25.      The factorization has the form
  26.         A = U' * U ,  if UPLO = 'U', or
  27.         A = L  * L',  if UPLO = 'L',
  28.      where U is an upper triangular matrix and L is lower triangular.
  29.  
  30.      This is the unblocked version of the algorithm, calling Level 2 BLAS.
  31.  
  32.  
  33. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  34.      UPLO    (input) CHARACTER*1
  35.              Specifies whether the upper or lower triangular part of the
  36.              symmetric matrix A is stored.  = 'U':  Upper triangular
  37.              = 'L':  Lower triangular
  38.  
  39.      N       (input) INTEGER
  40.              The order of the matrix A.  N >= 0.
  41.  
  42.      A       (input/output) REAL array, dimension (LDA,N)
  43.              On entry, the symmetric matrix A.  If UPLO = 'U', the leading n
  44.              by n upper triangular part of A contains the upper triangular
  45.              part of the matrix A, and the strictly lower triangular part of A
  46.              is not referenced.  If UPLO = 'L', the leading n by n lower
  47.              triangular part of A contains the lower triangular part of the
  48.              matrix A, and the strictly upper triangular part of A is not
  49.              referenced.
  50.  
  51.              On exit, if INFO = 0, the factor U or L from the Cholesky
  52.              factorization A = U'*U  or A = L*L'.
  53.  
  54.      LDA     (input) INTEGER
  55.              The leading dimension of the array A.  LDA >= max(1,N).
  56.  
  57.      INFO    (output) INTEGER
  58.              = 0: successful exit
  59.              < 0: if INFO = -k, the k-th argument had an illegal value
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSPPPPOOOOTTTTFFFF2222((((3333FFFF))))                                                          SSSSPPPPOOOOTTTTFFFF2222((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              > 0: if INFO = k, the leading minor of order k is not positive
  75.              definite, and the factorization could not be completed.
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.